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Under two different scenarios for the light scalar mesons, wc investigate the transition 
form factors of B{Bs) mesons decay into a scalar meson in the perturbative QCD approach. 
In the large recoiling region, the form factors arc dominated by the short-distance dynamics 
and can be calculated using perturbation theory. We adopt the dipole parametrization to 
' recast the dependence of the form factors. Since the decay constants defined by the scalar 

current are large, our predictions on the B ^ S form factors are much larger than the 
B P transitions, especially in the second scenario. Contributions from various light-cone 
distribution amplitudes (LCDAs) are elaborated and we find that the twist-3 LCDAs provide 
more than a half contributions to the form factors. The two terms of the twist-2 LCDAs give 
destructive contributions in the first scenario while they give constructive contributions in 
the second scenario. With the form factors, we also predict the decay width and branching 
ratios of the semileptonic B SlT> and B Sl^l~ decays. The branching ratios oi B ^ SlD 
' channels are found to have the order of 10^"' while those of B Sl^l~ have the order of 

10 ^. These predictions can be tested by the future experiments. 

^: 
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Although a number of scalar states have been discovered since long time ago, the underlying 
structure of scalar mesons has not been well established (for a review, see (H-|i|). In order to 
'sj" . uncover the inner structures, many different descriptions have been proposed such as qq, qqqq, 

I meson-meson bound states or even supplemented with a scalar glueball. It is very likely that they 

are not made of one simple component but are the superpositions of these contents. The different 
scenarios tend to give very different predictions on the production and decay of the scalar mesons 
OO ', which are helpful to determine the dominant component. Although intensive study has been given 

' to the decay property of the scalar mesons, the production of these mesons can provide a different 

■ unique insight to the mysterious structure of these mesons, especially their production in B decays. 

^ , In B meson decays, the energy release is much large and many channels involving a scalar meson 

H I in the final state are open. Since the first observation of the scalar meson /o(980) in three-body 

B meson decays B~ K~ fo{980) — > K~ {7r~^7r~) 0], the two collaborations, BaBar and Belle, 
have reported many studies on decays involving a scalar meson in the final state: the branching 
ratios and/or direct CP asymmetries are measured or set an upper limit [H]- Since much more 
interesting channels are still not observed at present, it is just the beginning of scalar meson study 
in B factories. Meanwhile, it is also necessary to provide more theoretical studies which are useful 
for future experiments. 

Theoretically, the studies on hadronic B decays are usually polluted by the nonperturbative 
QCD effect and predictions on the observables always suffer large uncertainties. Since there is 
only one hadron in the final state in semileptonic B ^ S decays, they receive less theoretical 
uncertainties. In these channels, the most challenging part in the calculation is the matrix element 
of the to scalar meson transition. In the region of small recoil, where is large, the form 
factors are dominated by the soft dynamics, which is out of control of perturbative QCD. However, 
in the large-recoil region where q'^ — > 0, roughly 5 GeV of energy is released. About half of this 
energy is taken by the light scalar meson, which suggests that large momentum is transferred in 
this process and the interaction is mainly dominated by the short-distance dynamics. Therefore 
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the perturbative QCD approach (PQCD) [gj] is expected to be applicable to B to scalar meson 
transitions in the large-recoil region. With the results obtained in the restricted region, one can 
extrapolate these form factors to the whole kinematic region by adopting some parametrization 
form for the form factors. 

This paper is organized as following: The distribution amplitudes and decay constants of the 
mesons are given in Section [III In Section IIIII we listed the formulae about the form factors and 
semileptonic decays. Section HV] are the discussion of the numerical results. The Appendix lAl lists 
out the useful functions for PQCD approach. 



II. CONVENTIONS AND INPUTS 

We will work in the rest frame of the B meson and use the light-cone coordinates. In the heavy 
quark limit the mass difference of b quark and B meson is negligible: mi, ruB- The masses of 
scalar mesons are very small compared with the b quark mass, we keep them up to the first order. 
Since the scalar meson in the final state moves very fast in the large-recoil region, we define the 
momentum of the scalar meson on the plus direction in the light-cone coordinates. The momentum 
of B meson and scalar mesons can be denoted as 

Pb = ^(1,1,0x), P5 = ^(r/, 0,0^). (1) 

Then for momentum q = Pb — Ps, there exists rj = 1 — q^/m^. The momentum of the light 
antiquark in B meson and the quark in scalar mesons are denoted as ki and k2 respectively(see 
FigllD: 

h = (0, ^xi, ki_L) , k2 = C^X2'n, 0, k2±) . (2) 

In the course of the PQCD calculations, the light-cone wave functions of the mesons are required. 
The B meson is a heavy-light system, and its light cone matrix element can be decomposed as 0] 



(2vr) 



^.)(^l) + ^^<^Sw(^l) 



(3) 

Pa 



where n = (1, 0, O-p) and v = (0, 1, Ot) are light-like unit vectors. There are two Lorentz structures 
in B meson light-cone distribution amplitudes, and they obey the normalization conditions: 



= , / T^Ms, iki) = 0, (4) 



with as the decay constant of B(^g^ meson. In principle, both the (pB^^^iki) and (pB^^^iki) 

contribute in B meson transitions. However, the contribution of 0^(3) (^i) is usually neglected, 
because its contribution is numerically small 0]. So we will only keep the term with (pB^^-^iki) in 
equation ([3]). In the momentum space the light cone matrix of B meson can be expressed as: 

^B(s) = ^^^^(s) +"^i?(.))75 (fci). (5) 



3 



Usually the hard part is independent of k'^ or/and k~ , so we integrate one of them out from 
(A;"'", k_L). With h as the conjugate space coordinate of k_L, we can express (^^f^j (2;, k^) in 
b-space by 

where x is the momentum fraction of the light quark in B meson. In this paper, we use the following 
expression for (pBf^sf{x,b): 



mn X 1 



(7) 



with Nsf^^s^ the normalization factor, which is determined by equation In recent years, a 

lot of studies for and decays have been performed by PQCD approach 0]. With the 
rich experimental data, the uJb in ([7]) is fixed as 0.40GeV. In our calculation, we adopt Uf, = 
(0.40 ± 0.05)GeV and fe = (0.19 ± 0.025GeV) for B mesons. For Bs meson, taking the SU(3) 
breaking effects into consideration, we adopt u>i)^ = (0.50 ± 0.05)GeV and /b^ = 0.23 it 0.03GeV[9j]. 

In the spectroscopy study, many scalar states have been discovered. Among them, the scalar 
mesons below 1 GeV, including /o(600)(c7), /o(980), K^{800){k) and ao(980), are usually viewed to 
form an SU(3) nonet; while scalar mesons around 1.5 GeV, including /o(1370), /o(1500)//o(1700), 
i^Q(1430) and ao(1450), form another nonet. There are two different scenarios to describe these 
mesons in the quark model. The first one(called scenario 1 in this paper) is the naive 2-quark 
model: the nonet mesons below 1 GeV are treated as the lowest lying states, and the ones near 
1.5 GeV are the first excited state. In this scenario, the flavor wave functions of the light scalar 
mesons are 

a = -^{uu + dd), fo = ss, 
v2 

Qq = ud, Oq = —^(uu + dd) , Qq = du, 
v2 

k"*" = us, = ds, = sd, K~ = su. (8) 

Here it's supposed that the a and fn (980) has the ideal mixing. However, the data of J/^l' decays 
doesn't favor /o(980) as a pure ss state[l0t], and it seems that a and /o(980) have a mixing like 

|/o(980)) = \ss) cos 9 +\nn) sin 6*, 

\a) = —\ss) sin 6 + \nn) cos 9, (9) 

with \nn) = -^{uu + dd) and 9 as the mixing angle. The above description has encountered 

several severe difficulties. For example, if the qq states have the quantum numbers J^'~^ = O"*"^, the 
corresponding masses are expected larger than that of the vector mesons. Studies on the mixing 
angle of a and fn (980) [ill] show that 9 tends to be not a unique value, which indicates that a and 
/o(980) may not be purely qq states. Based on these facts, the second scenario is proposed, where 
the nonet mesons near 1.5GeV are viewed as the lowest lying states, while the mesons below 1 GeV 
may be viewed as four-quark bound states. Because of the difficulty when dealing with four-quark 
states, we only do the calculation about the heavier nonet in this scenario. 



The decay constants of scalar mesons are defined by|10l] 



(5(^)1^27/^91 1 0) = fsPiM , {S\q2qi\0) = msfs- 



(10) 



4 



Because of the charge conjugate invariance, neutral scalar mesons cannot be produced by the vector 
current and thus 



f<^ = ffo = fa'^ 



0. 



(11) 



For other scalar mesons, the vector decay constant fs and scalar decay constant /s' (listed in Table 
|I]and[Tl]) is related by equations of motion Hsfs = fs, with fig = {Jy^mx (^) • mass of 

the scalar meson, and mi, m2 are the running current quark masses. Inputs of the scalar mesons 
in our calculation, include the decay constants, running quark masses in this paragraph and the 
Gegenbauer moments in the following, quote from [lo| . 

The definition of twist-2 light-cone distribution amplitude(LCDA) ^s{^) and twist-3 LCDAs 
<I>|(x) and for the scalar mesons can be combined into a single matrix element id ]: 



{S{Ps)\q{%q{zm 



-1 

-1 



with the normalization conditions 

'1 



dxe''^P-'{^S^s{x) + ms^Ux) + msiA A " 1)^5(2^)^7,(12) 



dx(f)s{x) 



fs 







/ dx(l)'s{x) = I 
Jo Jo 



1 

dx(l)g{x) 



fs 



(13) 



The LCDAs can be expanded in Gegenbauer polynomials as the following form: 

fs 



(t>s{x) 



^six) 
<fs{x) 



Is 
fs 



:6x(l — X) 
:6x{l — X) 



m=l 
^ oo 

- + Y^BMC'^H2x-l 



m=l 



^ ^ m=l 



d rs{x) 

dx 6 



Is d 
2^/2Wcdx 



.(l-x)[l + 5;^6„(/.)C7^/2(2x-l)] 



m=l 



(14) 
(15) 
(16) 



where BmifJ-), am{lj) and bmifJ-) are the Gegenbauer moments and clrJ"^^ and Cni^ are the Gegen- 
bauer polynomials. The values of Bm(n) are listed in Table HI and HH And the values of bmfJ. and 
amifJ') in scenario 2 are worked out in [20|, which is listed in Table Hill However, in the calculation 
in scenario 1, the asymptotic form of twist-3 LCDAs is used. 
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III. S FORM FACTORS AND SEMILEPTONIC DECAYS IN THE PQCD 

APPROACH 

A. A Brief Review of pQCD Approach 

The basic idea of pQCD approach is including the intrinsic transverse momenta of valence 
quarks in the calculation of the hadronic matrix elements. The transition matrix element (see Fig. 
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TABLE I: Decay constants fs (in unit of MeV) and Gegcnbauer moments at scale ^ = IGeV in scenario 1. 





fs 


Bi 


B3 


ao(980) 


365 ± 20 


-0.93 ±0.10 


0.14 ±0.08 


ao(1450) 


-280 ± 30 


0.89 ±0.20 


-1.38 ±0.18 


/o(980) 


370 ± 20 


-0.78 ±0.08 


0.02 ±0.07 


/o(1500) 


-255 ± 30 


0.80 ±0.40 


-1.32 ±0.14 


k(800) 


340 ± 20 


-0.92 ±0.11 


0.15 ±0.09 


if* (1430) 


-300 ± 30 


0.58 ±0.07 


-1.20 ±0.08 



TABLE IL Decay constants fs (in unit of MeV) and Gegcnbauer moments at scale /i = IGeV in scenario 2. 





Is 


Bi 


B3 


ao(1450) 


460 ± 50 


-0.58±0.12 


-0.49 ±0.15 


/o(1500) 


490 ± 50 


-0.48±0.11 


-0.37 ±0.20 


iff!; (1430) 


445 ± 50 


-0.57±0.13 


-0.42 ±0.22 



[T]) of B meson to a scalar meson(gi^2 component is supposed) can be expressed as the convolution 
of the wave functions ^b, find the hard scattering kernel Th, integrated over the longitudinal 
and transverse momenta of the valence quarks: 



M (X dxidx2 / 77-^77-^^>B(a;i,A-ix,PB,t)Tf/(xi,X2,A-i_L, A:2±,t)%(x2, A:2±,R,i)- (17) 
Jo J-00 (2vr)^ (27r)^ 

It's convenient to calculate the transition amplitude in coordinate space. Through the Fourier 
transformation, the above equation becomes 



pi poo 

M dxidx2 d\d'^b2'^Bixi,bi,pB,t)THixi,X2,bi,b2,t)<^s{x2,b2,Pi,t). (18) 

Jo J-00 

In principle, loop corrections to scattering kernel Th can be taken into consideration, which usually 
bring two types of infrared divergences in individual diagrams: soft and collinear. Soft divergence 
is generated when all the components of a loop momentum I go to zero: 

^^ = (/+,r,rT) = (A,A,A), (19) 



TABLE III; Gegcnbauer moments for the twist-3 LCDAs of scalar mesons at the scale /i = IGeV in scenario 
2(2^. 



state 


ai(xlO-2) 


02 




04 


fei(xl0-2) 


b2 


hi 


ao(1450) 





-0.33- 


-0.18 


-0.11 - 0.39 





- 0.058 


0.070 - 0.20 


iff!; (1430) 


1.8 - 4.2 


-0.33 - - 


-0.025 




3.7-5.5 


- 0.15 




/o(1500) 





-0.33 - 


0.18 


0.28 ~ 0.79 





-0.15 - -0.08 


8 0.044- 0.16 
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FIG. 1: Contributions to the form factors in the PQCD approach, where the cross denotes the weak vertex. 



with P expressed in the hght-cone coordinate. The cohinear divergence arise from the region where 
the gluon momentum is parallel to the massless quark momentum: 



(20) 



In both cases, the loop integration correspond to f d'^l/l'^ ~ log A, thus logarithmic divergences 
are generated. In perturbation theory, it has been shown order by order that these divergences 
can be separated from the hard kernel and obsorbed into meson wave functions using eikonal 
approximation 13] • When the soft and collinear momenta overlap, one also encounter double 



logarithm divergences, which can be resummed into the Sudakov factor and its expression is given 
in Appendix [Al 

The loop corrections to the weak decay vertex will generate another type of double logarithm. 
For example, the amplitude of the left diagram of Fig. [1] is proportional to When 
X2 0, additional collinear divergences are associated with the internal quark. The integration of 
the amplitude will produce double logarithm oAn^ X2, and the resummation of this type of double 
logarithm gives rise to Sudakov factor 5'j(x2)[13|], which is usually called jet function. The similar 
jet function St{xi) is generated after the resummation of the same type of double logarithm of the 
right diagram in Fig. [TJ The jet function decreases faster than any power of x as x — > 0, thus it 
kills the endpoint singularity effectively. The jet function has been parametrized in a form which 
is independent of the decay channels, twists and flavors [14 1. 

With the Sudakov factors included, the factorization formula of the form factor matrix element 
in pQCD approach is given by 



(21) 



M oc dxidx2 I d%id%2^B{xi,hi,pB,t)TH{xi,X2,bi,b2,t) 

Jo J-oo 

x^s{x2,b2,pi,t)St{xi)ex.p[-SB{t) - S2(t)]. 



B. Form Factors in the PQCD approach 



The form factors for B^^-^ — > S transition are defined by 



Ks{S{Ps)\q7f,75b\B^s){PB^^^)) = -i 



'^s{S{Ps)\qcrf,ub\B^s){PB^^^)) = -ie^uaisPiq' 



iPB,,,+Ps),- 



m 



mt 



mB^,) + ms 



Ks{S{Ps)\qcrf^ul5b\B(^s^{PB^^^)) = [q^,Psu- Ps^lqv 



2FT{q^ 



Fiiq') + 



ml 
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with q = Pb^^^-^ — Ps- ^5 is the flavor factor for the transition: \/2 for the component of uu in 
the state, for the component of dd in the state, 1 for the other states. In the 

large-recoil region, a hard gluon is required to kick the soft spectator antiquark to a fast-moving 
antiquark. Therefore, in this kinematics region, the form factors can be calculated perturbatively. 
The lowest order diagrams for the B(^g^ — > S transition are shown in Fig(TJ Carrying out the 
calculation under pQCD approach, we obtain the analytic formulae of the form factors nearby the 
q2 = 0: 

rl poo 

Poiv) = SvrCirm^ / dxidx2 / bidbib2db2(f>B{xi,bi) 
Jo Jo 

X { [r/(x2?? -r]- l)(l)s{x2) - r2T]{2x2 - l)(l)s{x2) + r2{2x2'n - 3t] + 2)(/)|(x2)] 
xhe{xi, (1 - X2)'i],bi, b2)as(tl)exp[-Sabitl)]Stix2) 

+2r2#|(x2)/ie(l - X2,xir], 62, bi)as{tl)exp[-Sab(tl)]St{xi)} , (25) 

Piiv) = ^nCprn^B I dxidx2 / bidbib2db2<i>B{xi-,bi) 
Jo Jo 

X { [(x2?? -1]- 1)05(X2) + r2(-22;2 + 3 - 2/vi)fs{x2) - ^2(1 - 2X2)(/)|(X2)] 
xheixi, (1 - X2)'i],bi, b2)as(tl)exp[-Sabitl)]Stix2) 

+2r2(ps{x2)he{l - X2,xir], 62, bi)as{tl)exp[-Sabitl)]St{xi)} , (26) 

Ft{7]) = 87rC_Fm|(l + r2) / dxidx2 / bidhib2db2(t)B{xi,bi) 

Jo Jo 

X { [r2(x2 - l)(t>s{x2) - (t)s{x2) + r2{x2 - 1 - 2/r])(f>^{x2)] 
xhe{xi, (1 - X2)'i],bi,b2)as(tl)exp[-Sabitl)]St{x2) 

+2r2(l)Ux2)he{l - X2,xir],b2,bi)as{tl)eM-Sab{tl)]St{xi)} . (27) 

With these formulae we calculate the form factors nearby = 0. Through fitting the results 
among the region < q^ < lOGeV^, we extrapolate them with the pole model parametrization 

^'^'^'^ " l-a{qyml) + b{qymiy' ^^^^ 
with a, b are the constants to be determined from the fitting procedure. 

C. Semileptonic Sj^j Meson decays 

The effective Hamiltonian for b uWi transition is 

n,ff{b ^ uWi) = ^Vubuj^il - 75)W7'^(l - (29) 
With the Hamiltonian, the q"^ dependant decay width ^ can be expressed as 



2 



dq^ ~ 192^3^1 (g2)2 Y g2 y 4^2 "^5 

(mf + 2g2)(g2 _ (^^ _ msf){q^ - {jub + Tasf)Fl{q^) + 3m2(m| - valfFl(ii^) [^0) 



with mi as the mass of the lepton. 
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The calculation of 6 ^ sl~^l~ transition is a bit complicated, because both the short-distance 
and long-distance contribution should be taken into consideration. The weak effective Hamiltonian 
is 

10 

'H^E = --^VMY.C,[^i)0,{^J), (31) 

^ i=l 

with the doubly CKM suppressed terms omitted. Ci{fi) are the Wilson coefficients and the local 
operators Oidj.) are given by [l^ 



Oi = 


{SaCa)v-A{cpbf3)v-A, O2 


= (Sa 


C/3)F-A(c/3&o)y-A, 


03 = 


{Saba)v^A y^(g/3g/3)y-A, 
<? 


04 = 


(gafc/3)y-A y^(g/3ga)y-A, 

Q 


05 = 


(gafca)y-A y^(g/3g/3)y+A, 


06 = 


{Sabi3)v-A'^{qpqa)v+A, 


07 = 






q 


09 = 


^(/7mO(^-7^(1-75)6), 


do = 





(32) 

where {qiq2)v-Aiq3qi)v-A = {qil''{'^-l5)q2){q37f,{'^-l)qi), and {qiq2)v-A{q3qi)v+A = (^i7^(l- 
75)92) (^37^(1 + 7)94)- In equation (f32l) . the term suppressed by in O7 is neglected. 
The amplitude for b — > sl^l~ transition can be decomposed as 

A(h ^ sl+l-) = 9l^^viVt,\clff{ii)[s^,{l - i,)b][iri] + do[5-7M(l - 75)6][/V75/] 

-2nn,C',ff{p) [sia^,^{l + 75)^'f7''^] } , (33) 

where s = jvr?^ and rhh = uib/mB, with as the b quark mass in the MS scheme. The 
long-distance and short-distance contributions are absorbed into the C^^-^ {n) and Cg'^-^(/i), with 

cfifi) = d(/.) + cu,(/^), 

Cf (/x) = Cgifl) + ypert(5) + rLD(s). (34) 

Fpert represents the perturbative contributions, and Ild is the long-distance part. The ^pert is 
given bvjl7l| 



Ypertis) = h{m„s)Co-h{l,s){iC3 + iCi + 3C5 + Ce) 

-h{0, s){C3 + 3d) + ^(3d + d + 3d + d), (35) 

with Co = d + 3d + 3d + d + 3d + d- The Wilson coefficients, listed in table HYI are 
given in the leading logarithmic accuracy. The long-distance part Ild, involving the contributions 
of — > SV{cc) resonances where V{cc) are the charmonium states, is neglected in this paper 
because of the lack of the experimental data. The corrections of the nonfactorizable effects of the 
charm quark loop to the 6 — > 57 transition at = are also neglected. And the absorptive part 
of 6 ^ 57 with neglecting the small contribution from Vtf,V^* is represented by the C^^^^ part in 

Cj^"^ , which is given by (for a complete expression of C^"^^ (fJ-), see [l^) 

C^^,^(/i) = ias[^V^^/^HGi{xt) - 0.1687) - 0.03d(^*)], (36) 



9 



TABLE IV: The values of Wilson coefficients Ci{mi,) in the leading logarithmic approximation in Standard 
Model, with mw = 80.4GeV, mt = 173.8GeV, rub = 4.8GeV.[i3 



Ci 


C2 C3 


C4 C5 




Cj Cg Cio 


1.119 


-0.270 0.013 


-0.027 0.009 - 


-0.033 


-0.322 4.344 -4.669 



TABLE V: Form factors for _B ^ S* in scenario 1. The errors arise from the mrcertainties of hadronic 
parameters of _B(s) meson(/f, and cuh), Aqcd, scalcs(ig) and the Gegenbauer moments of scalar mesons. 



Fo(0) = Fi(0) Ft{0) a{Fo) b{Fo) a{F,) b{F,) aiFr) HFt) 



B 

B ao(980) 
B k(800) 



0.39 



^-0.06 
+0.10 
'0.08 



27+"-°^ 

-0.06 



0.29 
0.45 
0.29 



0.07 
-0.06 
+0.11 
0.10 
+0.07 
0.07 



'-'•"0-0.07 

r, 79+0.08 

'^-0.03 

71 +0-04 

' ^-0.08 



-0 11+000 

"•^^-0.13 



-0.16: 
-0.12: 



-0.12 
-0.00 
-0.02 
-0.12 



1.61 
1.68 
1.65 



+0.04 
-0.06 
+0.03 
0.06 
+0.06 
0.04 



0.56+o°4 



0.62 
0.59 



-0.10 
+0.01 
0.10 

+0.08 
0.04 



1.67+0-05 



1.70 

1.69 



0.05 
+0.06 
0.03 
+0.06 
0.05 



n 62+° °^ 
'J-'J^-o.oo 

63+° " 

"-'•O'5-O.Ol 

+0.08 
0.06 



0.65] 



B ^ /o(1370) 
B ao(1450) 
B A'o*(1430) 



-0.30 
-0.31 
-0.34 



+0.08 
-0.09 
+0.08 
0.09 
+0.07 
0.09 



-0.39tE;;}? 
-0.41+°i° 



-0.44 



"^-0.12 
+0.10 
-0.11 



+0.07 
-0.02 
+0.13 
-0.02 



0.70 

070+0.13 

+0.04 
0.04 



0.72^ 



-0.24 
_0.26+"-24 

5+0.04 
0.05 



+0.15 
0.05 
0.24 
-0.00 



1 RQ+0.09 

^•"'J-0.05 
^_g3+0.08 



-0.18: 



1.65 



-0.04 
+0.04 
-0.07 



^■^•5-0.08 

0.53+°-^3 



0.57 



-0.06 
+0.08 
-0.14 



1 60+°°^ 
^•"'J-0.04 

1.62+°04 



1.61 



-0.07 
+0.04 
-0.05 



n 50+° °^ 

0.54+0.03 



0.52 



-0.13 
+0.05 
-0.06 



BO ^ /o(980) 
B", k(800) 
50 ^ /o(1500) 
A* (1430) 



r, or+0.09 

"••J "3-0. 07 

0.29+0.07 



-0.26 
-0.32 



0.06 
f0.09 
-0.08 
+0.06 
-0.07 



-0.10 
-0.08 
+0.07 
■^-0.06 

_0.34+o.io 



0.40 
0.31 



-0.41 



-0.10 

+0.08 
-0.09 



0.73 
0.66 
0.72 
0.69 



0.04 
-0.06 
0.07 
-0.03 
+0.14 
0.08 
+0.05 
-0.03 



10+0.02 

-0 17+0-11 
^•^ ' -0.00 

,0.20+010 



-0.21 



-0.10 
+0.11 

-0.03 



1 66+°°*^ 

^^•""-0.05 
+0.03 
-0.05 
+0.13 
-0.03 
+0.06 
-0.03 



1.62 
1.61 
1.62 



^•""-0.05 

n 5fi+o oo 

'-'■OO-0.09 
0.48+0.27 



1 70+0 06 

J^- '"-0.04 



1.68] 



0.52 



-0.02 

+0.14 
-0.04 



0.05 
-0.03 

1.60+0.06 



1.62 



-0.04 
+0.01 
-0.06 



^•""-0.05 

62+° i° 
'^•"J^-0.01 

0.48+0.09 



0.56 



-0.04 
+0.00 
-0.16 



with Gi{xt) = 8(k-i)^ + 4(4-1)4 ' ^ = oisimw)/as{fJ.) and xt = m4/ml^. 
The dependant width oi B ^ Sl~^l~ is given by 



1024m|7r5 



q^ 



2 2 y iriB + ms 



^F,{q') 



-Fi{q')^ 



l^^esmtFTiq^) 



'7 



ms + ms 



+ 



with A = (m^ — q 



4A 

2 



m 



ruB + ms 



+ 4 



m;Cio(m| - m|)Fo(g^) 



(37) 



2 \2 



2 „2 



Arrigq 



IV. NUMERICAL RESULTS AND DISCUSSION 



Form Factors 



Our results of the B ^ S form factors are listed in table |V] and IVIi The errors for the form 
factors in those two tables arise from the uncertainties of hadronic parameters of B^^^-^ meson(/B and 
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TABLE VI: Form factors for S ^ 5 in scenario 2, with the same error sources as the data in Table fVl 



Fo(0) = Fi(0) Ft(0) aiFo) b{Fo) a(Fi) 6(Fi) aiFr) biFr) 



B ^ /o(1370) 
B ao(1450) 
B K^{USO) 



0.63l°-23 
0.68+°i^ 



0.60 



0.15 
+0.18 
0.15 



0.76t°;?^ 

j+0.25 
-0.19 



^- '"-0.11 
^■"^-0.08 



-0 14+002 
-0.21+°°6 



0.78"' 



0.68 



+0.07 
-0.05 



-0.18 



-0.02 
+0.06 
-0.01 



1 60+015 

^•"'J-0.05 
1 70+0.12 
^- ''-'-0.07 

'"-0.07 



"•^•5-0. 09 

0_70+o.i6 



0.65 



-0.11 
+0.10 
-0.10 



cq+0.07 
""3-0.05 
RS+006 

+0.07 
0.04 



68"* 



n f^7+0 07 
"■^'-0.07 

"■"-'^-0.02 
+0.11 
0.02 



o.6r 



> /o(1500) 



0.60 
0.56 



+0.20 
0.12 

+0.16 
-0.13 



0.82 



+0.30 
-0.16 
,+0.22 



0.65 



+0.04 
0.10 
7+0.06 



'-'•'^-0.17 "•"'-0.07 



-0.22 

-0.17 



+0.07 
-0.02 
+0.01 
-0.07 



1.76 
1.69 



+0.13 
-0.08 
+0.08 
-0.07 



0.71 
0.63 



+0.20 

-0.08 
+0.09 
-0.10 



1.71 



+0.04 r, aa+O.OG 

-0.07 "-""-O.IO 

1 fis+0 06 n R-i+om 

^■0^-0.06 U-O;i_o.08 



^b)^ AQCD(0.20GeV-0.30GeV), factorization scales(see Eqs. ()Aip ) and the Gegenbauer moments of 
scalar mesons. A number of remarks will be given in order. 

• Compared with transitions of B meson to pseudoscalar mesons, vector mesons and axial- 
vector mesons our predictions on i? ^ S" form factors in scenario 2 are obviously 
larger, which is caused mainly by the large decay constants(/5) of the scalar mesons. For 
example, the form factor i<b(0) of B meson to pion transition is about 0.230] with O.lSlGeV 
as the decay constant of pion, while the B meson to ao(980) transition in scenario 1 has 0.39 
as its corresponding form factor, whose decay constant is more than two times larger than 
pion. 

• In Table IVl the form factors of — > cr are smaller than those of i? ^ ao(980). Because the 
same decay constant and Gegenbauer moments for these two particles are used in the cal- 
culation, the differences are caused by the mass differences between ao(980) and cr(0.98GeV 
for ao(980) and 0.513GeV for cj[l|). In scenario 1, there are small differences between k(800) 
and /o(600) in masses(0.672GeV for k(800)), decay constants and Gegenbauer moments. 
Besides, the contribution from twist-2 LCDA of k(800), which is proportional to fs , is 
too small to give sizable differences. Thus the B ^ a and B — > k(800) have nearly the 
same form factors as shown in Table |Vl Comparing the form factors of -B — > ft:(800) with 
B^ — > k(800) in Table IVl one can find that the differences between B and B^ mesons affect 
little. Therefore, the large differences between the form factors of B^ — > k(800) and those of 
B — > ao(980) are mainly due to the large difference between the scalar meson masses. 

• The form factors of B to heavier nonet transition in scenario 1 are negative, while the others 
are positive. The reason is that the decay constants(/5) of the heavier nonet in scenario 1 
have opposite signs to the others, which is clearly shown in Table HI and HIl 

• As we can see from the table |V] and IVIl the predictions in scenario 2 are larger than the 
corresponding ones in scenario 1 roughly by a factor of 2 in magnitude. In order to show 
how these large differences are generated, we take the form factor Fq{Q) as an example and 
list contributions from different terms in LCDAs in Table IVTiy Data is given with asymptotic 
forms of twist-3 LCDAs are adopted in both scenario 1 and scenario 2, because the terms 
with Gegenbauer moments bring so small effects, which is discussed in the following, that 
they can't change the argument). The contributions from the two twist-3 LCDAs 4>'^g and (frg 
are given in the first two columns. The numbers in the column denotes the contributions 
from the Gegenbauer moments Bi in twist-2 LCDAs. It is also similar for the fourth i?3 
column. The last column collects the total contributions to the form factors. The different 
inputs between in scenario 1 and in scenario 2 are the decay constants and Gegenbauer 
moments. If only twist-3 LCDAs are taken into account, the form factors will be proportional 
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to the decay constant. Since the decay constants fs in S2 are (typicahy 60%) larger than 
those in SI in magnitude, the form factors are accordingly larger. The (pg term give much 
larger contributions than the (fg term. Contributions from the Gegenbauer moments of the 
twist-2 LCDAs sizably enhance the form factors in S2 but not too much in SI. For B to 
scalar meson transitions in scenario 1, the Bi terms provide contributions with the same sign 
with the twist-3 terms, while the terms with B^ have the opposite sign. Thus the two terms 
of the twist-2 LCDAs give destructive contributions to the total form factors in SI. The 
situation is different in S2, although the two Gegenbauer moments are small in magnitude, 
they give constructive contributions and induce much larger form factors. 

We also investigate the contributions from terms with Gegenbauer moments in twist- 
3 LCDAs, and find that the effects brought by these moments are not large. Taking 
B — > /o(1370) transition as an example, a comparison between the cases with and with- 
out these contributions is given in Table IVIIIi We can see that most of the results are 
changed by less than 10%. 



Compared with our previous study on i? — > /o, i^'g (1430) transitions [2l|, |2^, the predictions 
for the form factors given in the present work are a bit smaller. The main reason is that 
different values for the threshold resummation parameters c have been used. Moreover, the 



form factors in this paper are larger than those obtained in other approaches or models [27l. 
28l . 2^ . 30]. As a result, the branching ratios of the semileptonic decays are larger, which is 
discussed in the following. 

As we have mentioned in the introduction section, the experimentalists have already provided 
many investigations on nonleptonic B decays involving a scalar meson in the final state. Among 
these decays, the so-called color-allowed tree-dominated processes can be directly utilized to esti- 
mate the B ^ S form factors, under the hypothesis of factorization. For example, the B^ — > aQiT~ 
decay amplitude in the factorization scheme is expressed as: 



AiB^^a+TT~) = -^m|/^Fo^^"°{KbKrd[ai+«4 + aio-r^(a6 + a8)] 

+^cfeKd["4 + flio - r^iae + as)]} , (38) 
where a, is the combination of Wilson coefficient 

ai = C2 + Ci/3, a2 = Ci + C2/3, 
ai = Ci + Ci+i/N, (i = 3,5,7,9), 

a^ = Ci + Ci^i/N, {i = 4, 6, 8, 10). (39) 

ai 1, and it has small uncertainties. Although there are large uncertainties for 03-010, the 
combination of Wilson coefficients satisfies: 

ci » maxJa3„io]. (40) 

If only the branching ratios are concerned, contributions from the penguin operators (03-10 terms) 
can be safely neglected and thus 

^(^o^a+vr-) = ^ml/.Fo^^'^oKfcKVi- (41) 

If the partial decay widths are well determined experimentally, these results will directly constrain 
the B to scalar meson transition form factors. The upper bounds for B — > aovr are given as(in unit 
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TABLE VII: Contributions from different LCDAs to the B ^ S form factor Fq in scenario 1 (SI) or 
scenario 2 (S2). The contributions from the two twist-3 LCDAs 0^ and cj)^ are given in the first two 
columns. The numbers in the column' Bi denotes the contributions from the Gegenbauer moments Bi in 
twist-2 LCDAs. It is also similar for the fourth column. The last column collects the total contributions 
to the form factors (Data is given with asymptotic forms of twist-3 LCDAs adopted in both scenario 1 and 
scenario 2). 





<j>'s (j)^ Bi Bs Total 


B ao(1450) 


SI 
S2 


-0.21 -0.05 0.14 -0.19 -0.31 
0.35 0.08 0.15 0.11 0.69 


B iq{U30) 


SI 
S2 


-0.22 -0.05 0.10 -0.18 -0.34 
0.33 0.07 0.14 0.09 0.62 


5° /o(1500) 


SI 
S2 


-0.17 -0.04 0.11 -0.16 -0.26 
0.32 0.08 0.13 0.09 0.61 


S," ^ A- (1430) 


SI 
S2 


-0.19 -0.05 0.09 -0.17 -0.32 
0.27 0.07 0.14 0.09 0.58 



TABLE VIII: Form factors for B /o(1370). The first line and the second line are the results with and 
without contributions from the terms with Gegenbauer moments in twist-3 LCDAs respectively. 



Fo(0)-Fi(0) 


Ft{0) 




HFo) 


a(Fi) 


6(Fi) 


aiFr) 


HFt) 


67+"- 1^ 




0.701^1 
n 71+0 02 

L"- ' -L-0.07 


-0 14+002 
-0 12+°°° 


1 60+°-i5 
1 64+°°'' 


fS3+°-i^ 

'-'•O'J-0.09 

fS7+° °'' 


1 fi^+°°5 


0.571°:°? 



of 10-6): 

Bn{B 0^(980)77+) < 3.1, 

BTZiB 0^(1450)^+) < 2.3, (42) 

where the daughter BF has taken to be 100%. Since the scalar mesons ao(980) and ao(1450) have 
vanishing decay constants in the isospin limit, the branching ratios of — > 7r+ are very small 
and one expects the relation: B1Z{B — > a^7r+) = BTZ{B^ a^T^'). Compared with the branching 
ratio of B^ 77+77" (in unit of 10"^) 

57^(S ^ 77+77") = 5.16 ±0.22, (43) 

results provide the upper bound for the B ^ form factors: 

Fo{B ao(980)) < 0.78Fo{B ^ 77) = 0.18, Fo{B ao(1450)) < 0.67Fo(5 ^ 77) = 0.15, (44) 

where as an rough estimation, we have taken Fq(B 77) = 0.233. Compared with our results 
in Table |V] and IVIl one can see our results have exceeded the present experimental upper bound. 
Despite of that, it does not mean our predictions are ruled out by the data, as the daughter decay 
is not taken into account in the derivation for the experimental bound. Our predictions will be 
confronted with the real bound in the future, whenever the daughter decay of oq is well studied. 
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TABLE IX: The total branching ratios for the b — > ulfi in scenario l(Unit:10 ''). The errors are estimated 
with errors from the form factors. 





B Seve{ ^iv^i) B StVt 


B- ^ a 
B° a+(980) 
B- ^ /o(1370) 
B° ^ a^{U50) 


o.8il°J? o.silS:?^ 

1.84lJ-°^ l.OltnJ 

0.291°:!^ o.is;^^ 

O.GTl^i^ 0.28tg;l^ 


Bs -> K+(800) 
Bs ^ /Vo*^(1430) 


1.42l°:i o.88lS:g 



TABLE X: Same as Table |TX] except m scenario z. 





B- 


^ /o(1370)ei/e( //i/^) 




/o(1370)Ti>, 


This work 




1 r;c + 1.53 
-■^•^0-0.65 




67+0-68 






a([(1450)ez?e( 






This work 




Q oc;+2.36 
'-'•^^-1.36 




1 Q9+0.97 
-'^•'-'^-0.57 


LCSR[27] 








0.63l°:i 




5s - 


^ K*+{UiQ)eVe{ny^,) Bs- 


> /-i:o + (1430)Ti>^ 


This work 




^■^3-1. 05 




L09l°i? 


LCSR[2a 




^•«5-0.4 




0.521°:?^ 


QCDSR[2a 




U.OO_Q 24 







B. Decay widths and branching fractions 

With the form factors at hand, one can directly obtain the partial decay width through Eq. (j30|) 
and Eq. pzp . Since masses of electrons and muons are very small compared with in most 
kinematic region of the semileptonic decays, they will not produce large effects and are neglected 
in this work. In Fig. [2]and[3l we give our predictions on the partial decay width of -6(5) — > Sl~i>i 
{I = e,/i) and -B^^) ~^ St~i>i^ respectively. The diagrams in Fig. [Hand Fig. [5] are similar but for 
the -6(5) Sl'^l' {I = e, ^) and B^^g^ — > St^t~ decays. In Fig. dj there exists a small discontinuity 
in each diagram, which is caused by the discontinuities in functions h{rhc,s) and h{l,s) in Eqs. 
(jSSp . When Z = r in Fig. [5l the discontinuities in the diagrams disappear, because the origins of 
q'^ axes become 4m^ which is large enough to ensure that the variation of does not pass the 
discontinuities in the h{mc,s) and h{l,s) functions. 

The results for the total branching ratios are collected in Table IIXI [Xl IXII and IXIIl with the 
errors estimated with the errors of the form factors. One can find that the branching ratios with 
r lepton(s) in the final state are smaller than the ones without r lepton(s), because the large mass 

of T lepton(s) makes the phase space much smaller. In Table IIXI is smaller than 

two when the scalar meson belongs to the light nonet. While for the heavy nonet mesons, the 
value of this ratio is larger than two. The reason is that more energy is released when the final 
state is a light meson, and thus the effect of nir on the phase space is not so evident. In Table 1X1 
and IXIIl we also list the predictions in light-cone sum rules(LCSR) and QCD sum rules(QCDSR), 
which are smaller than our predictions. The reason is that we have bigger form factors. Taking 
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FIG. 2: Partial decay widths of the semileptonic B Slv decays as functions of . Diagram a-d denote 
the B- (cr,aJ(980),/o(1370),aJ(1450))r 'vi in scenario 1 respectively; Diagram e-f denote the B 
(/o(1370), a(]'(1450))Z~i/i in scenario 2 respectively; Diagram g: Bs k+(800)/~z?; in scenario 1; Diagram 
h: Bs — > K^'^ {1A30)1~ Di in scenario 1; Diagram i: Bg K^^ {1430)1^ Di in scenario 2. 
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FIG. 3: Partial decay widths of the semileptonic B Stv decays as functions of . Diagram a: 
<7T~Di in scenario 1; Diagram b: _B" {980)t'' Di in scenario 1; Diagram c: B^ /o(1370)t~/?; in 

scenario 1; Diagram d: i?° aJ(1450)T~P; in scenario 1; Diagram c: B~ fo{1370)T~ i>i in scenario 2; 
Diagram f: B^ — > {14:50)t~ Di in scenario 2; Diagram g: Bg — > k^{800)t~Di in scenario 1; Diagram h: 
Bs /\q^(1430)t~i>; in scenario 1; Diagram i: Bg A'Q^(1430)T~t'i in scenario 2. 
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7.5 10 12.5 

(Diagram f) 



FIG. 4: Partial decay widths of the semileptonic B Sl^l {I — e, fi) decays as functions of q^. Diagram 
a: B^ K^l^l^ in scenario 1; Diagram b: B^ K^^ {1430)1'^ l~ in scenario 1; Diagram c: B^ 
K*~{U30)l+l- in sccnario_2; Diagram d: i?° /o(980)/+/- in scenario 1; Diagram e: B° /o(1500)/+/- 
in scenario 1; Diagram f: /o(1500)/"*'/- in scenario 2; 



TABLE XI: The total branching ratios for the b — > sl^l in scenario l(Unit:10"'^) with the same error 
sources as Table HXl and 1x1 





B - 


5*6+6 ( /X+/i ) 


B -> St+t- 


B- K- 
B- i^*-(1430) 




o 1 q + 1.73 
0.10_]^ 21 


'-'■OO-0.25 

2.00tJ;7? X 10-2 


5° ^ /0(980) 
5° /o°(1500) 




c oi +3.23 
"J-^1^-2.06 

1- '^-0.94 


2.2ll};^? X 10-2 



B'^ — > a^(1450)e fg as an example, the form factors that contribute are i^oC?^) and Fi{q'^), with the 
relationship Fo(0) = Fi(0). K)(0) for 5° aj"(1450) in scenario 2 in this paper is O.GS^j};}^, while 
the corresponding value in [27| is 0.52 it 0.10. As a rough estimation, supposing that corresponding 
form factors in these two papers have analogical evolution with respect to q^, the branching ratio 
in this paper should be (0.69/0.52)^ ~ 1.7 times larger. 
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FIG. 5: Partial decay widths of the semileptonic B St^t^ decays as functions of . Diagram a: B^ ^ 
K^T^T^ in scenario 1; Diagram b: Kq~ {1A3Q)t^t^ in scenario 1; Diagram c: B^ A'Q^(1430)r+r- 

in scenario 2; Diagram d: B'^ /o(980)r+r- in scenario 1; Diagram c: 5° /o(1500)r+r^ in scenario 1; 
Diagram f: 5° — > /o(1500)r+r- in scenario 2; 



TABLE Xll: Same as Table IXll except in scenario 2. 





B- A'*"(1430)e+e-( 


-) B- /i*-(1430)T+r- 


This work 


9.78llig 


6.29t^-^^ X 10-2 


LCSR[27] 


c 7+3.4 
0- ( _2.4 


9.8t^2.4 ^ 


LFQM[29] 


1.63 


2.86 X 10-2 


QCDSR[30] 


2.09-2.68 


(1.70- 2.20) X 10-2 




BO^/0(1500)e+e-(A*+M") 


i?o ^ /o"(1500)r+r- 


This work 




n 13+"i2 

^•^'-'-0.06 


LCSRf27] 


c 0+2.3 


n i2+o-°® 



V. CONCLUSIONS 

In this work, we have studied the B ^ S form factors in the PQCD approach under two different 
scenarios for the scalar mesons. In scenario 1, both of the light and heavy nonet are described as the 
qq state while in scenario 2, we have only studied the heavy nonet. Due to the large decay constant 
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fs, we have found that most of our predictions are larger than those for the B ^ P transition form 
factors, especially in scenario 2. Contributions from various LCDAs are explicitly specified. Due to 
the large masses of oo(1450), Kq (1430), /o(1500), their twist-3 LCDAs have provided more than one 
half contributions to the form factors in both scenarios. In scenario 1, the two Gegenbauer moments 
Bi,B3 for the twist-2 LCDAs have different signs and they give destructive contributions to the 
form factors; while in scenario 2, although the two Gegenbauer moments are small in magnitudes, 
they give constructive contributions and induce larger form factors. Contributions from terms with 
Gegenbauer moments in the twist-3 LCDAs are also investigated, and we find that these terms 
do not give large changes. We also study the semileptonic B Slv and B Sl~^l~ decays, 
including the partial decay width and the integrated branching fractions. Branching ratios of the 
semileptonic B Slv decays are found to have the order of 10~^, while branching fractions of 
the B Sl'^l' decays have the order of 10~^. Compared with results in the previous studies, our 
predictions are a bit larger which is caused by larger form factors. These predictions will be tested 
by the future experiments. 
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APPENDIX A: PQCD FUNCTIONS 

In this part, we collect the functions which are essential in the PQCD calculation. 



t\ = max(tc\/ (1 - X2)rimB, 1/fei, I/62), = lnax{tc^/xrfjmB, l/h, I/62), (Al) 
with tc = 1 for the calculation of the central values and tc = 0.75-1.25 for error estimation. 



/le(xi,X2, 61,62) = Ko{y/xiX2mBbl) 



9{bi - b2)Ko{y/x^mBbi)Ioiy/x^mBb2) 



+0(62 - bi)KQ{,/x2mBb2)h{^/x2mBbi) 



(A2) 



s.(x)^ ™/!V^ [-a-)r^ 

V7rr(l + c) 

with c = 0.4. The Sudakov factor in Eqs. (l25p -Eqs. ()27p is given by 

Sab{t)=SB{t)+Ss{t), 

where 

SB{t) = s(^x,^M)+ll^^^^l,{as{fl)), 

Ss{t) = .(x2^,62)+.((l-X2)^,62)+2^*^ ^7,K(/i)), 



(A3) 
(A4) 

(A5) 
(A6) 
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with the quark anomalous dimension 7^ = —as/iT- The exphcit form for the function s{Q,b) is: 



s{Q,b) 



2/3i 
4/33 



oln - 
' b 



2/3i 



b + 



ln(2g) + 1 ln(25) + 1 



+ 



b 

8/33 



4/32 4/3i 
ln2(2^) - ln2(26) 



In 



where the variables are defined by 

q^ln[Q/{V2A)], S^ln[l/(6A)], 
and the coefficients A^*) and Pi are 



/3: 



33 - 2n 



^(2) 



12 

67 7r2 10 



153- 19n/ 
24 



3 27"/+3A'"<2'™>- 



(A7) 



(A8) 



(A9) 



is the number of the quark flavors and is the Euler constant. We will use the one-loop 
running coupling constant, i.e. we pick up only the four terms in the first line of the expression for 
the function s{Q,b). 
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